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I. INTRODUCTION 
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. The nucleon energy momentum tensor (EMT) form factors [1] contain valuable information on the nucleon structure. 
' They carry information on, e.g., how the quark and gluon degrees of freedom share the total momentum, and angular 
O momentum of the nucleon [2], or on the distribution of strong forces inside the nucleon [3]. The first is known 
from deeply inelastic lepton nucleon scattering experiments. The latter two can be deduced from generalized parton 
[ distribution functions (GPDs) [4-7] accessible in hard exclusive reactions [8-13], see [14-18] for reviews. 

Nucleon EMT form factors were studied in lattice QCD calculations [20-24]. In principle, lattice QCD provides a 
[ rigorous and model-independent approach to compute the nucleon EMT form factors. In practice, however, present 
. day technics and computing power allow to simulate on lattices "worlds" with pions of typically TOtt ^ 400 MeV. The 
situation is extected to improve in the future, see [25] for status reports on selected topics. 

For the time being, however, it is necessary to use chiral extrapolation in order to relate lattice results to the 
' real world situation. Chiral perturbation theory (xPT) provides a model-independent tool for that, and the chiral 
, behaviour of the nucleon EMT form factors was studied in [26-28] . Experience with chiral extrapolations of other 
C — ■ nucleon properties indicates that xPT is applicable up to the lowest presently available lattice values of m-^ [29-31] 
] although the issue is not yet settled [32]. 

. In this situation it is worth looking on what one can learn about the chiral behaviour of nucleon properties from 
rS i' other effective approaches, e.g. the "finite range regulator" (FRR) approach. There chiral loops are regulated by 
' suitably chosen vertex form factors in order to phenomenologically simulate the effects of the pion cloud which has a 
finite range due to ^ [33, 34]. However, model calculations [35-37] are equally of interest in this context. 
, A phenomenologically successful and theoretically consistent model for the description of nucleon properties at the 
■ physical point and in the chiral limit, is the chiral quark soliton model (CQSM) [38, 39]. This model describes the 
nucleon as a soliton of a static background pion field in the limit of a large number of colours Nc, and hence provides 
a particular realization of the general large- TVc picture of the nucleon [40]. The CQSM describes numerous nucleonic 
?H ' properties without adjustable parameters - including among others form factors [41-43], usual parton distribution 
functions [44-48] and GPDs [49-54] — within an accuracy of (10 — 30)% as far as those quantities are known. 

That it is possible to extend the CQSM to the description of the nucleon at large pion masses was shown in [37], 
where the model was demonstrated to provide a good description of lattice data on the m7r-dependence of the nucleon 
mass Mjv up to m,r ^ 1-5 GeV. An important prerequisite for that is that the CQSM formally contains the correct 
heavy quark limit result for the nucleon mass Mjy [37]. 

In this work we present a study of the m^-dependence of the nucleon EMT form factors in the CQSM up to pion 
masses as large as 1 GeV. The present study extends the study of Ref. [55] where the nucleon EMT form factors were 
studied at the physical point and in the chiral limit, and its purpose is threefold. 

First, we provide an important supplement for the study in Ref. [37]. There soliton solutions were obtained 
numerically for model parameters corresponding to pion masses in the range < < 1.5 GeV. Here we provide a 
cross check demonstrating that the numerical solutions found in [37] correspond, in fact, to stable solitons. 

Second, with the results obtained for large we are in the position to confront the model predictions for the 
nucleon EMT form factors directly to lattice QCD results [20-24]. In view of the early stage of art of the experimental 
situation of hard exclusive reactions, such a comparison provides the only presently available test for our results. 

Third, though the model — as discussed in detail in [37] — cannot be used as a quantitative guideline for the chiral 
extrapolation, our study still allows to gain several interesting qualitative insights with this respect. 

The note is organized as follows. In Sec. II we introduce the nucleon EMT form factors and discuss their properties. 
In Sec. HI we briefly review how the nucleon EMT form factors are described in the CQSM. In Sees. IV and V 
we describe respectively the model results for the densities associated with the form factors and the form factors 
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themselves. In Sec. VI we compare the model results with lattice QCD data, and in Sec. VII we discuss which 
qualitative observations from our study could be of interest in the context of the chiral extrapolation of lattice data. 
Sec. VIII contains the conclusions. A remark on different notations for the EMT form factors is posed in App. A. 



II. FORM FACTORS OF THE ENERGY-MOMENTUM TENSOR 

The nucleon matrix element of the symmetric EMT of QCD is characterized by three scalar form factors [1]. The 
quark and gluon parts, and T^, of the EMT are separately gauge-invariant and can be parameterized as [2, 3], 
see App. A for an alternative notation, 



^ ^ ' Mn 2Mjv 



+ ^ '^'^''^■^ 



u 



ip)- (1) 



Here the nucleon states and spinors are normalized by {p'\p) = 2p°(27r)^(5(^^(p' — p) and u{p)u{p) — 2Mn, and spin 
indices are suppressed for brevity. The kinematical variables are defined as P = {p + p')/2, A = {p' — p), t = A^. 
The form factor c{t) accounts for non-conservation of the separate quark and gluon parts of the EMT, and enters the 
quark and gluon parts with opposite signs such that the total (quark+gluon) EMT is conserved. 

The nucleon form factors of the EMT are related to the second Mellin moments of the unpolarized GPDs H^{x, ^, t) 
and E-f{x,^,t) as (we use the notation of Ref. [16]) 

1 1 

[dxx^Hf{x,i,t)=M^{t) + U?{t)e , fdxxJ2iH^'+Ef){x,^,t)=2jQ{t). (2) 
-if -if 

where ^ denotes the so-called skewedness parameter [2] . The sum rules in Eqs. (2) are special cases of the so-called 
polynomiality property of GPDs [14]. The second sum rule in (2) provides the possibility to access J'^(O), i.e. the 
total (spin-horbital angular momentum) contribution of quarks to the nucleon spin, through the extraction of GPDs 
from hard exclusive processes and extrapolation to the unphysical point t = 0. The sensitivity of different observables 
to the total quark angular momenta was investigated in model studies [16, 56]. For gluons there are analog definitions 
and expressions. Suffice to remark that the full GPDs contain far more information [19]. 

The form factors of the EMT in Eq. (1) can be interpreted [3] in analogy to the electromagnetic form factors [57] in 
the Breit frame characterized by A*^ = 0. In this frame one can define the static energy- momentum tensor for quarks 

T^Ar, = ^ / ^ exp(zAr) {p' , S'\f^M\P, S) , (3) 

and analogously for gluons. The initial and final polarization vectors of the nucleon S and S' are defined such that 
in the respective rest-frame they are equal to (0, s) with the unit vector s denoting the quantization axis for the spin. 

The components of T^(r, s) and e''''^rjT^(r, s) correspond respectively to the distribution of quark momentum 
and quark angular momentum inside the nucleon. The components of {T^ — ^SikT^ ){r,s) characterize the spatial 
distribution of "shear forces" experienced by quarks inside the nucleon. The respective form factors are related to 

jQ^t) + ^jQ'{t) = ld're-^e^^'s,rjTg{r,s), (4) 
d?{t) + jd?\t) + '^drit) = /d^re-^r^T^^W (rV-^^-) , (5) 

where the prime denotes derivative with respect to the Mandelstam variable t. Note that for a spin-1/2 particle only 
the T°''-components are sensitive to the polarization vector. Note also that Eq. (6) holds for the sum Tqo = + Tgg 

with M2{t) = M^it) + M2^{t) and J{t) and di{t) defined analogously, but not for the separate quark and gluon 
contributions - since otherwise the form factor c{t) would not cancel out. 
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The form factor M2{t) at f = is connected to the fractions of the nucleon momentum carried respectively by 
quarks and gluons. More precisely 



1 1 



(7) 



where /f (a;) = H°'{x,0,0) are the unpolarized parton distributions accessible in inclusive deeply inelastic scattering. 

The form factors Af^''^(<), J'^''^'{t) and df''^{t) arc rcnormalization scale dependent (the indication of the renor- 
malization scale n is suppressed for brevity). Their quark+gluon sums, however, are scale independent form factors, 
which at f = satisfy the constraints. 



M2(0) = ^ yd3rroo(r,s) = l, 
J(0) = Jd\s'^''sirjTok{T,s) = ^, 



di{0) = I d^r T,,(r) (rV^ - j 6^^^ ^ d, , (8) 

which mean that in the rest frame the total energy of the nucleon is equal to its mass, and that its spin is 1/2. The 
value of di is not known a priori and must be determined experimentally. However, being a conserved quantity it 
is to be considered on the same footing as other basic nucleon properties like mass, anomalous magnetic moment, 
etc. Remarkably, di determines the behaviour of the D-term [58] (and thus the unpolarized GPDs) in the asymptotic 
limit of rcnormalization scale ^ ^ oo [16]. 

The form factor di{t) is connected to the distribution of pressure and shear forces experienced by the partons in 
the nucleon [3] which becomes apparent by recalling that Tjj(r) is the static stress tensor which (for spin and 1/2 
particles) can be decomposed as 

^ii(r) = sir) - ^ S,,>j + p{r) S,j . (9) 

Hereby p{r) describes the radial distribution of the "pressure" inside the hadron, while s(r) is related to the distribution 
of the "shear forces" [3]. Both are related due to the conservation of the EMT by the differential equation 

2 5£(r) ^ MO ^ W ^ _ (10) 

3 dr r dr 

Another important consequence of the conservation of the EMT is the so-called stability condition 



oc 

/■ 



dr r^p{r) = . (11) 

Let us review briefly what is known about di — which in terms of the pressure or shear forces is given by 

di = -^MN j d^r s(r) = ^ Mn j d^r p{r) . (12) 

For the pion it can be calculated exactly using soft pion theorems with the result df^ = —M^^ [58], see also [59]. Also 

for the nucleon < was found in calculations in CQSM [49, 55, 60]. For the nucleon the large- A^c limit predicts the 
flavour-dependence -|- | » jd" — | [16]. Lattice calculations [20-24] confirm this flavour dependence and yield 
a negative rf^, see Sec. VL In a simple "liquid drop" model di can be related to the surface tension of the "liquid" 
and comes out negative [3]. Such a model is in particular applicable to nuclei and predictions from this picture [3] 
were confirmed in calculations assuming realistic nuclear models [67]. In Ref. [55] it was conjectured on the basis of 
plausible physical arguments that the negative sign of d\ is dictated by stability criteria. This conclusion, however, 
remains to be proven for the general case. 
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III. NUCLEON EMT FORM FACTORS IN THE CQSM 



The effective theory underlying the CQSM was derived from the instanton model of the QCD vacuum [61-63] which 
assumes that the basic properties of the QCD vacuum are dominated by a strongly interacting but dilute instanton 
medium, see the reviews [64]. In this medium light quarks acquire a dynamical ("constituent") quark mass due to 
interactions with instantons. At low momenta below a scale set by w 600 MeV, where pav denotes the average 
instanton size, the dynamics of the effective quark degrees of freedom is given by the partition function [65, 66] 



(13) 



Here we restrict ourselves to two flavours and neglect isospin breaking effects with m = ?n„ = denoting the current 
quark mass, while U = exp{iT'^TT'^) denotes the chiral pion field with U^^ = exp{ij5T'^n°'). The dynamical mass is 
strictly speaking momentum dependent, i.e. M = M{p). However, in practical calculations it is convenient to work 
with a constant M = Ad{Q) ~ 350 MeV following from the instanton vacuum [64], and to regularize the effective 
theory by means of an explicit regularization with a cutoff of 0{p~^) whose precise value is fixed to reproduce the 
physical value of the pion decay constant f-^ = 93 MeV. In this work we will use the proper-time regidarization. 

The quark degrees of freedom of the effective theory (13) correspond to QCD quark degrees of freedom up to 
corrections which are small in the instanton packing fraction Pav/Rav ^ where i?av denotes the average separation 
of instantons. The same parameter suppresses the contribution of gluon degrees of freedom [63]. 

The CQSM is an application of the effective theory (13) to the description of baryons [38, 39]. While the Gaussian 
path integral over fermion fields in (13) can be solved exactly, the path integral over pion field configurations can be 
solved only in the large- iVc limit by means of the saddle-point approximation (in the Euclidean formulation of the 
theory). In the leading order of the large- A^c limit the pion field is static, and one can determine the spectrum of the 
one-particle Hamiltonian of the effective theory (13) 

H\n) = En\n) , H ^ -i-f^-f'^dk + -j^MlP' + 7"m . (14) 

The spectrum of (14) consists of an upper and a lower Dirac continuum, distorted by the pion field as compared to 
continua of the free Dirac-Hamiltonian (given by H with replaced by 1), and of a discrete bound state level 
of energy -Eicv, if the pion field is strong enough. By occupying the discrete level and the lower continuum states each 
by Ac quarks in an anti-symmetric colour state, one obtains a state with unity baryon number. The soliton energy 



E^i[U] = Nc 



Elev + ^ {En — Eno ) 
En<0 



(15) 



is a functional of the pion field. It is logarithmically divergent, see e.g. [42] for explicit expressions in the proper- 
time regularization. By minimizing Eso\[U] one obtains the self-consistent solitonic pion field Uc- This procedure is 
performed for symmetry reasons in the so-called hedgehog ansatz 7r"(x) = e" Pcir) with the radial (soliton profile) 
function Pc{r) and r = |x|, e,. = x/r. The nucleon mass is given by Esoi[Uc\. 

In the large- Ac limit the path integral over U in Eq. (13) is solved by evaluating the expression at Uc and integrating 
over translational and rotational zero modes of the soliton solution. In order to include corrections in the 1/Ac- 
expansion one considers time dependent pion field fiuctuations around the solitonic solution. In practice hereby one 
restricts oneself to time dependent rotations of the soliton field in spin- and flavour-space which are slow because the 
corresponding soliton moment of inertia 



N, 



_£ \ " 



I) {I 



E — E 

-'-'777, -*-'r) 



(16) 

reg 



is large, / = 0{Nc). It is logarithmically divergent and has to be regularized. In (16) one has to sum over occupied 
("occ") states n which satisfy En < £^iev, and over non-occupied ("non") states m which satisfy Em > Eiev 

The model expressions for the nucleon EMT form factors in the effective theory (13) were derived explicitly in [55]. 
The gluon part of the EMT is zero in the effective theory (13), because there are no explicit gluon degrees of freedom. 
(So we omit the index Q when discussing the model results in this and in Sees. IV, V but restore it later.) 

Consequently, in the model the quark part of the EMT is conserved by itself, and the form-factor c{t) in Eq. (1) 
vanishes [55] . The model expressions for the other form factors read 
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^'^^^ ~ ^'^'^^ = aZ;^ /^''" P"^''^ Mr^t) (17) 
d,it) = i^/d3r,(r)M!:^ (18) 

J{t) = 3|d^p.(r)^ii^, (19) 

with the Bessel functions jn{z) = ^^j^ and ji{z) = -~jl^{z). The Fourier transforms of the form factors, which are 
radial functions and to which we refer as "densities" in the following, are defined as 



n,occ 



reg 



P{r) = ^^C(r)(7°7P)<^n(r) 



3 

n.occ 



= -^E^"'^->Kr)(2p'' + (i;„ + i;,)7V)<^n(r) ^"'^''^'^ 



n,occ \ / J • 



(20) 
(21) 

(22) 



j.non 



The expressions in (20, 21, 22) are logarithmically UV- divergent. Here we use the proper-time method to regularize 
them, see [55] for explicit expressions in this regularization. In Ref. [55] analytical proofs were given that 

• the stability condition (11) is satisfied in the model, 

• the form factors satisfy the constraints at t = in Eq. (8), and 

• the same expressions for EMT form factors follow in the model from unpolarized GPDs via the sum rules (2). 

Notice that the latter is a special case of the "polynomiality property" of GPDs [14] satisfied in the CQSM [51, 52]. 
The field-theoretic character of the model is a crucial prerequisite which allows to formulate and analytically prove 
such and other [44] general QCD requirements. This in turn provides important cross checks for the theoretical 

consistency of the approach. 

For the numerical calculation we employ the so-called Kahana-Ripka method [71], whose application to calculations 
of the nucleon EMT form factors in the CQSM is briefly described in Ref. [55]. The use of the proper- time regu- 
larization has the advantage (over, e.g., the Pauli-Villars method [72]) that it is possible to include explicitly chiral 
symmetry breaking effects due to a finite current quark mass m in the effective action (13). 

In Ref. [37] it was shown that it is possible to obtain soliton solutions and compute nucleon masses for current 
quark masses up to rn = C(700MeV) which corresponds to pion masses up to TOtt = 0(1.5 GeV). What provides a 
certain justification for the application of the model up such large m is the fact that the model formally contains the 
correct heavy quark limit result for the nucleon mass [37]. The proof that in the limit m — > mg, where mg is the 
heavy quark mass, the nucleon mass tends to Mjv NcTUq given in [37] is formal because in this proof it was taken 
for granted that stable soliton solutions do exist up to such large pion mass values. 

Therefore, it was of importance in Ref. [37] to demonstrate numerically the existence of soliton solutions for large 
up to, at least, = 0(1.5 GeV). These soliton solutions were found by using a standard iteration procedure for 
the calculation of the self-consistent profile function Pc{r) described in detail, e.g. in [42]. Here, as a byproduct of our 
study of EMT form factors, we will be in the position to provide an important and valuable cross check. Namely, do 
the pressures p(r) computed with the respective large-m„ soliton profiles really satisfy the stability condition (11)? 
The answer is yes, see below Sec. IV D. 

Before discussing the m^-dependence of the densities (20, 21, 22) and the form factors (17, 18, 19) we have to 
establish which model parameters are allowed to vary and which are kept fixed while is varied. Here we shall use 
the choice of Ref. [37] to keep M = 350 MeV and /^r = 93 MeV fixed. Then the proper-time cutoff is adjusted such 
that for a given m and one reproduces the physical value of (m and m,r are related to each other in the effective 
theory (13) by a relation which for small m corresponds to the Gell-Mann — Oakes — Renner relation, see [37]). 

This way of parameter handling in the model was found to provide a good description of lattice data on the variation 
of Mjv(m,r) with m,r, once one takes into account the generic overestimate of the nucleon mass in the soliton approach 
[69] . However, the above way of parameter handling is just one possible choice and other choices are possible as well. 
An investigation whether other choices of parameter handling yield equally satisfactory results will be presented 
elsewhere. 
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IV. THE DENSITIES OF THE EMT 



In this Section we shall focus our attention on the densities (20, 21, 22) which are interesting objects by themselves, 
before we discuss the form factors (17, 18, 19) in the next Section. The study of the densities will enable us to address 
the question whether the model provides a satisfactory description of the nucleon in (fictions) worlds with pion masses 
up to 1.2 GeV. As we shall see, the answer is yes. 



A. Energy density 



The energy density psir) is just T°'^(r) in the static EMT (3), and is normalized as / d'^r ^^(r, to^) — MN{mT,) for 
any m^r where we explicitly indicate the pion mass dependence. Mffim^) as function of was studied in [37]. 

Fig. la shows psir) as function of r for pion masses in the range < < 1.2 GeV. In the following we focus on 
the region > 140 MeV, and include only for completeness the results for m-^ < 140 MeV discussed in detail in [55]. 

In the physical situation with = 140 MeV the energy density in the center of the nucleon is 1.7GeVfm~ or 
3.0 X lO^^gcm"^. This corresponds roughly to 13 times the equilibrium density of nuclear matter. As ttitt increases 
Pe{0) becomes larger and reaches p_b(0) = 9.5GeVfm~^ = 17 x 10^^ g cm~'^ at m,r = 1.2 GeV. At the same time, 
with increasing nij^ the fall-off of psif) at large r becomes stronger, as can be seen in Fig. lb. 
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FIG. 1: (a.) The energy density Pe{^) of the nucleon as function of r for different pion masses, (b) The normalized energy density 
47rr^PE(r)/A/]v as function for different pion masses. The curves are normalized such that one obtains unity upon integration over r. 
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FIG. 2: The pion mass dependence of (a) the energy density Pe{0) in the center of the nucleon, and (b) the mean square radius of the 
energy density (r^) defined in (23). The crosses indicate the physical point. 
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These observations mean that with increasing m^r the nucleon becomes "smaher" . To quantify this statement we 
consider the mean square radius of the energy density defined as 



. l^f^ , (23) 
J d'^r pE{r) 



which decreases with increasing TOtt- The pion mass dependence of P£;(0) is shown in Fig. 2a, see also Table I 
where many results are summarized. We observe an approximately linear growth of psiO) with mj. The pion mass 
dependence of {r'^) is shown in Fig. 2b, see also Table I. Up to m^r ^ 400 MeV we observe a roughly linear decrease 
of (r|;) with m,r which proceeds at a slower rate for TO,r ^ 400 MeV. (Throughout we choose a linear or quadratic in 
rriTr presentation of the mjr-dependence of the quantities — depending on which one is more convenient.) 

The above observations can be intuitively understood. With increasing the range of the "pion cloud" decreases. 
This results in a less wide spread nucleon. The above observations are also consistent with what one expects from 
the heavy quark limit point of view. The heavier the constituents building up a hadron, the smaller is the size of that 
hadron. Thus, the model results for Pe{t) a-re in agreement with what one expects for increasing m^. 

We remark that, being a chiral model, the CQSM correctly describes the behaviour of (r|;) in the chiral limit [55]. 



B. Angular momentum density 

The angular momentum density pj{r) is related to the T'^^ components of the static EMT as p,j{r) — t^^^ SiXjT^ki'^)- 
Fig. 3a shows pj(r) as function of r for different pion masses. For any rriTr we find that p,j{r) oc at small r, it 
reaches then a maximum around r ~ (0.3 — 0.4) fm and goes slowly to zero at large r. 

As increases we observe that the density pj(r) becomes larger in the small r region at the price of decreasing in 
the region of larger r. The increase in one and decrease in another region of r (as TOjt is varied) occurs in a precisely 
balanced way, because pj{r) — in contrast to the energy density — is always normalized as /d^r pj(r, rnvr) = Jn = 5, 
independently of m-^ . These observations can be understood within the picture of a rigidly rotating soliton as follows. 
For large pion masses the "matter" inside the soliton is localized more towards its center, as we have observed above, 
such that the inner region of the soliton plays a more important role for its rotation. As m^r decreases, and hence the 
range of the pion cloud increases, the energy density in the soliton becomes more strongly delocalized, and then the 
"outer regions" play a more and more important role for the rotation of the soliton. 

These findings can be quantified by considering the mean square radius (rj) of the angular momentum density 
defined analogously to (23). Fig. 3b shows (j'j) as function of 171-^, and we see that (rj) decreases with increasing m-^. 
Notice that in the chiral limit pj{r) oc at large r such that (rj) diverges [55]. Our numerical results for (rj) in 
Fig. 3b indicate this effect. 




FIG. 3: (a) The angular momentum density pj{r) as function of r for different pion masses m^. (b) The mean square radius (r^) of the 
angular momentum density as function of m,r ■ 
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C. Pressure and shear forces 



Next we turn to the discussion of the distributions of pressure and shear forces, p{r) and s{r), which are related to 
the T**^ components of the static EMT. 

Figs. 4a and b show the distributions of pressure p{r) and shear forces s(r) as functions of r for different m^. For 
all m,r the distributions of pressure and shear forces exhibit the same qualitative behaviour. The pressure takes at 
r = its global maximum, decreases monotonically becoming zero at some point rg till reaching its global minimum at 
some point rp_ min, and decreases then monotonically tending to zero but remaining always negative. The distribution 
of shear forces is never negative. It starts at a zero value at r = 0, increases monotonically till reaching a global 
maximum at some point rs^max, and decreases then monotonically tending to zero. 

The positive sign of the pressure for r < ro corresponds to repulsion, while the negative sign in the region r > ro 
means attraction. This is intuitive because in the inner region we expect repulsion among quarks due to the Pauli 
principle, while the attraction in the outer region is an effect of the pion cloud which is responsible for binding the 
quarks to form a nucleon [55] . 

With increasing the pressure p{0) in the center of the nucleon increases. At the same time also the absolute 
value of its (negative) minimum increases. Also the maximum of s(r) becomes larger with increasing ttItt, while the 
characteristic positions rg, rp^min and rs_ max move towards smaller r, see Fig. 5 for m7r-dependence of p(0) and rg. 

These observations can be understood as follows. The repulsive forces in the center of the nucleon increase as a 
response to the higher density at larger ttt-tt. At the same time the size of the nucleon decreases requiring stronger 
binding forces — and a "movement" of characteristic length scales of p{r) and s(r) towards the center. At any m^r 
repulsive and attractive forces are precisely balanced due to (11), see next Sec. IVD. 
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FIG. 4: The distributions of (a) pressure p{r) and (b) the shear forces s{r) as functions of r for different pion masses. 
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FIG. 5: The pion mass dependence of the pressure p(0) in the center of the nucleon, and of the position ro at which p{r) vanishes. 
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FIG. 6: (a) r^p(r) as function of r for different pion masses. The integrals over the regions where r'^p(r) is respectively positive and 
negative cancel each other within numerical accuracy of better than 1%. This shows how the stability condition in (11) is satisfied, see 
text, (b) The pressure p{0) as function of the energy density Pe{0) in the center of the nucleon. Some pion masses to which the respective 
values refer are indicated. The physical point is marked by the cross. 



D. Stability 



While the densities, PEii") and pj{r), are normahzed with respect to Mjv and Jn, for the pressure p{r) the 
corresponding anafogon is the stabiUty criterion (11). In Ref. [55] it was proven analytically that (11) is satisfied in 
the model — provided one evaluates the pressure with the self-consistent profile, i.e. with that profile which for a 
given TO,r provides the true minimum of the soliton energy (15). 

For given model parameters the soliton profiles are obtained by means of an iteration procedure which is described in 
detail for example in [42] . The profiles used here were computed in Ref. [37] where a good convergence of the iteration 
procedure was observed. However, what precisely means that the convergence of the iteration was good? In other 
words, how to test the quality of the numerical results? One could, for example, slightly modify the obtained profiles 
and check that they yield larger soliton masses than the respective true self-consistent profile. But the probably most 
elegant method is provided by the stability criterion (11). If, and only if, we found the soliton profile which truely 
minimizes the soliton energy (11), the pressure computed with that profile will satisfy (11). 

One way to check to which numerical accuracy our results satisfy (11) is as follows. Let us consider r'^p(r) as 
function of r, see Fig. 6a, and compute the integrals from to ro and from ro to oo.^ We obtain 



r 2.614 MeV °f r -2.630 MeV \ dr r'^v{r)\ f 0-31 % for = 140 MeV 

dr r2p(r) <^ 3.737 MeV / dr r^p{r) = I -3.748 MeV i.e. ^^ -f^ = <^ 0.15% for = 500 MeV 

U.856MeV, l-3.861MeV, Jo dr r |p(r)| (o.07% for = 1200 MeV, 

(24) 

and see that the stability criterion is satisfied to within a satisfactory numerical accuracy. 

Finally, we may test another sort of stability. We may ask the question how do pressure and energy density in 
the center of nucleon depend on each other for varying mjr. In fact, with /0£;(r, TO,r) and p(r, m^r) at hand, we may 
eliminate at r = and express p(0) as function of pe{Q). This is shown in Fig. 6b which demonstrates how the 
center of the nucleon responds to changes of to^ . Understanding the center of the nucleon for a moment as a "medium 



which is subject to variations of the external parameter m7r" we observe that for any rriTr we have 



dpjPE 
dpE 



> 0. This is 



a criterion for stability of a system which must respond with an increase of pressure if the density is increased. 



^ The numerical calculations are carried carried out in a finite spherical volume — here of the size D = 12 fm. For most quantities the 
densities decay fast enough at large r such that it is suflicient to integrate up to r = D. This is what we did in Eq. (24). However, for 
certain quantities given by integrals over the densities weighted by a higher power of r the integrands may happen not to be negligibly 
small at large r in particular for small m„ < 140 MeV. Then it is necessary to explore the analytically known large-r asymptotics of 
the densities, see [55], and to include the contribution of the regions r > D not covered in the numerical calculation. Examples of such 
quantities are di or (r j). The latter is divergent in the chiral limit, see Sec. IV B. 
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V. RESULTS FOR THE FORM FACTORS 



Fig. 7 shows the form factors of the EMT as functfons of t for \t\ < 1 GeV^ for different to^. All EMT form factors 
(with the exception of J{t) and di{t) in the chiral limit, see below) can be well approximated by dipole fits of the kind 

It is instructive to compare within the model the EMT form factors to the electromagnetic form factors — for example 
to the electric form factor of the proton Gsit) [41]. Interestingly, J{t) and Gsit) show a similar t-dependence. But 
M2{t) falls off with increasing \t\ slower than Gsit), while di{t) exhibits a faster fall off, see [55] for more details. 

The dipole masses of the different form factors exhibit different rriTr-dependences, see Fig. 8a and Table I. For all form 
factors the dipole masses increase with increasing m^. It is an interesting observation that the dipole masses of M2{t) 
and J{t) exhibit for m^r ^ 140 MeV to a good approximation a linear dependence on m^r. But the mTr-dependence of 
the dipole mass of di{t) follows a different pattern. We shall comment more on that in Sec. VII 

That the dipole approximation for J{t) and di{t) fails in the chiral limit, is due to fact that the slopes of J{t) and 
di{t) dX t — diverge in this limit [55]. For J{t) this is clear because its derivative at t = is related to the mean 
square radius of the angular momentum density as J'(0) = and (rj) diverges for — > 0, see Sect. IV B. 

The slope of di (t) at t = becomes infinitely steep in the chiral limit because it is related as d'l (0) = /d^r r^p{r) 
to the pressure which behaves as p{r) cx ^ at large r in the chiral limit. More precisely, d[(t) (x 1/V— t at small t 
in the chiral limit. For small but non-zero the derivative di(0) exists and is proportional to 1/m^. These results 
hold both in the CQSM [55] and in xPT [27]. The derivative of the form factor M2{t) at t = is finite for any TOtt. 

M2{t) and J{t) are normalized at t = for any m^r as M2(0) = 2J(0) = 1 [55]. In the CQSM these constraints are 
consistent for they mean that entire momentum and spin of the nucleon are carried by quark degrees of freedom. The 
numerical results satisfy these constraints within a numerical accuracy of better than 1 %, see Figs. 7a and 7b. 

In contrast, no principle fixes the normalization of the form factor di{t) at i = neither in the model nor in QCD. 
For all rn^ we find di = di(0) < 0. This condition has been conjectured to be dictated by stability requirements [55]. 
Fig. 8b shows the m,r-dependence of di which is rather strong. This is due to the fact that di receives a large leading 
non-analytic contribution proportional to ttIt^. (The "non-analyticity" refers to the current quark mass m oc mj.) 
The chiral expansion of di reads 

di(m^) = di -f fc + . . . (26) 

64 tt/^ 

o 

where di denotes the chiral limit value of c?i, and the dots indicate subleading terms in the chiral limit. Since the 
limits Nc ^ oo and — > do not commute [78, 79] one has in Eq. (26) k = 1 for finite [27] and fc = 3 in the 
large- iVc limit [55]. The latter corresponds to the situation in the CQSM. 

It is interesting to observe that the leading non-analytic term in the chiral expansion of di in Eq. (26) dominates 
the chiral behaviour of di up to the physical point, see Fig. 8b. But for larger TO,r higher orders in the chiral expansion 
become important, and change the qualitative m7r-behaviour of di(m^). We shall come back to this point in Sec. VII. 




FIG. 7: The form factors of the energy momentum tensor M2{t), J(t) and d\{t) as functions of t for different pion masses. In the chiral 
hmit J{t) and di (t) exhibit infinitely steep slopes at t = 0, and di (0) takes the value —3.46 which does not fit on the scale in Fig. 7c. 
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FIG. 8: (a) The dipole masses as defined in (25) for the form factors M2{t), di{t) and J(t) from the CQSM vs. m,r- Notice that di(t) and 
J{t) cannot be approximated by dipoles for m,r = 0. (b) The constant di as function of m,r at low scale. Solid line: The exact result from 
the CQSM. Dotted line: The chiral expansion of (ii(m^) to linear order in m,r according to (26). The square marks the physical point. 



Finally, we discuss mTr-dependence of the mean square radius (r|n) of the 
trace of the total EMT operator given due to the trace anomaly [80] by 



f/ = ^ F^'^F^, + (1 + 7„0 J2 ma^a^a 



(27) 



Let F{t) denote the form factor of the operator (27) which is normalized as 
F(0) = 1. Its slope at t = defines 



(4) = (4) 



p/ = 6F'(0) which can be related as 
12 di 



(28) 



to (r^) and di, see [55]. Fig. 9 shows how (rp) depends on the pion mass. 
In the chiral limit (rp) is the mean square radius of the operator F^'^F^y 
and its large value there is in contrast to what is known about the mean 
square radii of other gluonic operators [81]. 



<rp> infm^ 




0.2 0.4 0.6 0.8 in GeV 



FIG. 9; The mean square radius {rp) of the 
trace of the total EMT operator vs. mn. 



TABLE I: The pion mass dependence of different quantities computed in the CQSM using proper time regularization: the energy density 
in the center of the nucleon pe{0), the mean square radii (r^) and (r^) , the pressure p{0) in the center of the nucleon, the position rg 
of the zero of the pressure defined as p(ro) = 0, the constant di, the dipole masses of the form factors M2{t), J{t) and di{t), and the 
mean square radius (rp) of the trace of the EMT operator. In the chiral limit J(t) and di{t) have infinitely steep slopes at t = 0, see 
text. In these cases dipole fits do not provide useful approximations and are undefined (labelled by "undef." in the Table). The results 
for m.^ < 140 MeV from Ref. [55] are included for completeness. 





Pe(0) 


(-1) 




p(0) 


ro 


di 


dipole 


masses .A/^ip 


in GeV for 




McV 


GoV/fm^ 


fm^ 


fm^ 


GcV/fm^ 


fm 




Ahit) 


J{t) 




fm= 





1.54 


0.82 


oo 


0.195 


0.59 


-3.46 


0.867 


undef. 


undef. 


1.04 


50 


1.57 


0.76 


1.88 


0.202 


0.59 


-3.01 


0.873 


0.692 


0.519 


0.95 


140 


1.70 


0.67 


1.55 


0.232 


0.57 


-2.35 


0.906 


0.745 


0.646 


0.81 


300 


2.14 


0.53 


1.11 


0.298 


0.54 


-1.81 


0.990 


0.844 


0.872 


0.62 


500 


3.10 
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0.77 
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0.51 


-1.66 
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0.46 
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0.32 


0.59 
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0.49 
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1.214 


0.37 
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6.86 


0.26 


0.48 


0.553 


0.46 


-1.55 


1.334 


1.237 


1.337 


0.29 


1200 


9.53 


0.22 


0.38 


0.597 


0.42 


-1.47 


1.473 


1.390 


1.492 


0.24 
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VI. COMPARISON TO LATTICE RESULTS 

It is instructive to compare the results for the form factors of the EMT to lattice QCD data [20-24] . Presently, this 
offers actually the only available test for the model results. For the comparison it is necessary to evolve the model 
results from a low initial scale /Xq ~ 0.6 GeV to typically /x ~ 2GeV in the lattice calculations which we shall do to 
leading logarithmic accuracy. Under evolution the quark (flavour-singlet) form factors mix with the corresponding 
gluon form factors. We set the latter to zero at the initial scale. 

In this context it is worthwhile recalling that in early parameterizations of unpolarized parton distributions fi{x), 
the gluon (and sea quark) distribution(s) and thus M^' were assumed to be zero at a low initial scale [73] . One success 
of these approaches was that they were able to explain the observation « 0.5 at /i^ ^ 5 GcV^. I.e. starting with 
M2 = 1 and M2 = at a low scale, which is the situation in the CQSM, it is possible to reproduce the observation 
M2 ~ 0.5 at several GeV^. However, with the advent of more data (especially at low x) it became clear [74] that 
non-zero gluon (and see quark) distributions are required already at low initial scales. Modern parameterizations 
performed at low scales require a sizeable gluon distribution and M2 « 0.3 [75]. This is not in disagreement with 
the instanton picture where twist-2 gluon operators are suppressed with respect to quark operators by the instanton 
packing fraction which is numerically of order 30 % [63] . Thus in some sense the phenomenologically required "portion 
of gluons" is within the accuracy of the model [46]. With these remarks in mind we conclude that the CQSM result 
M2 = 1 at the low scale of the model is in agreement with phenomenology — within the accuracy of this model. 

Next, before we compare the results from the model to those from lattice QCD, let us confront the lattice results 
[20-24] with predictions from the large- A^c limit. In this limit one has for |t| ^ Mf^ independently of the scale [16] 

(A"+A'^)(t)» |(A"-A'')(i)| , |(B"~B'')(t)| > |(B" + B'')(t)| , |(C" + C"')(t)| » |(C;"-C'^)(t)| , (29) 

V ' ^ V ' ^ V ' ^ V ' V ' ^ V ' 

O{N0) 0{N-'') 0{N,) O(N0) 0{N^) 0{N,) 

cf. App. A for the explanation of the notation. For completeness let us quote that the gluon form factors satisfy 

A^{t) = 0{N^^), B^{t) = 0{N^), C^{t) = 0{Nl), (30) 

which is the same large- A^c behaviour as the corresponding quark-flavour-singlet form factors. 

Remarkably, although in the real world A^c = 3 does not seem to be large, nevertheless lattice data [20-24] reflect 
the large- A^c flavour dependence of the quark form factors (29). In fact, large- A^c relations of the type (29, 30) are 
observed to be satisfied in phenomenology [76] and serve within their range of applicability as useful guidelines [77]. 

The soliton approach is justified in the large- A^c limit [40] and satisfies general large- A^c relations of the type (29). 
The observation that the lattice results [20 24] are compatible with (29) is therefore an encouraging prereqiusite for 
our study. 

Let us first compare the model results to the lattice data computed by the LHPC and SESAM Collaborations [22]. 
There unquenched SESAM Wilson configurations on a 16^ x 32 lattice at /3 = 5.6 with k = 0.1560 were used. This 
corresponds to = (896 ± 6) MeV and a lattice spacing of a = 0.098 fm with physical units fixed by extrapolating 
the nucleon mass. The form factor A^+'^{t) = AQ{t) = M^{t) was computed omitting disconnected diagrams at a 
scale n = 2 GeV for < |f] < 3.1 GeV^. (The different notations are discussed in App. A.) The lattice data for A'^{t), 
which can be fit to the dipole form, are shown in Fig. 10a. In Fig. 10a we also show the CQSM results evolved to the; 
same scale for \t\ < 1 GeV^. We observe that the model results agree with the lattice data [22] to within 15% which 
is within the accuracy to which the CQSM results typically agree with phenomenology. 

Fig. 10b shows the form factor B^^'^if) = B'^ (t) = 2jQ{t) - M^{t) from Ref. [22] which was computed in the range 
—t G [0.59, 3.1] GeV and was found consistent with zero within the statistical accuracy of the simulation. Also in the 
CQSM we find B^{t) close to zero — in reasonable agreement with the lattice data, see Fig. 10b. Notice that in the 
model, at the low scale, we have _B^(0) = since (0) = 2J'^{0) (and equal unity), see Sees IVA and IVB. This 
implies a vanishing quark contribution to the " gravitomagnetic moment" of the nucleon [53] conjectured in [70]. The 
smallness of B^{t) implies that to a good approximation 2J^{t) w M^(t). We shall come back to this point below. 

Next, let us compare to the quenched lattice results by the QCDSF Collaboration [23] which were obtained using 
non-perturbatively 0{a) improved Wilson fermions on a 16'^ x 32 lattice at (3 = 6.0 using varying values of k. In 
Figs. 11a and lib and 11c the form factors Al+'^it) = A'^{t) = M^{t) and B'l^+'^{i) = B'^{t) = 2jQ{t) - M^(t) and 
C^+'^(t) = C"3(t) = idi(t) are shown for = 640 MeV [23]. The results refer to a scale of = 2 GeV, and cover 
the regions -t € [0, 2.8] GeV^ for A^it) and [0.6, 2.8] GeV^ for B^{t) and C"3(t). For comparison we plot the CQSM 
results at the corresponding scale and value of m^. Also here we observe a satisfactory agreement with the lattice 
data within model accuracy and/or within the statistical accuracy of the lattice data. 
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FIG. 10: The form factors of the energy momentum tensor A^+''{t) = AQ{t) = (t) and B^+'^it) = BQ{t) = 2jQ{t) - M^{t) as 
functions of t at the scale fi = 2 GeV in a world with pion mass m,.^ = 900 GcV. (See App. A for the discussion of different notations.) 
The lattice data points are by the LHPC and SESAM Collaborations [22]. The solid curves are the results from the CQSM obtained here. 
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FIG. 11: The form factors of the energy momentum tensor A^'^'^{t) = AQ{t) = (t) and Bl^+'^{t) = B^ (t) = 2.]Q{t) - M^{t) and 
Cl^+'^(t) = C'^{t) = ^df{t) as functions oft at the scale fi = 2GeV in a world with pion mass = 640 GeV. (For notations see App. A.) 
The lattice data points are from the QCDSF Collaboration [23]. The solid curves are the results from the CQSM obtained here. 



Notice that the presence of explicit kinematical factors of C(A) in the case of B'^{t), and of ©(A^) in the case of 
C^{t), see Eq. (Al) in App. A, practically amplifies the statistical errors of the form factors — as can be seen in 
Figs. 10 and 11. This is also the reason why no results for C'^{t) were presented in the exploratory study of Ref. [22]. 

The value of 2j'^{t) = {A'^ + B^){t) at t = is of particular interest since it gives the percentage of the total 
nucleon spin carried by quarks [5] . The lattice calculations yield 

0.60 ± 0.07 at /i = 1.76 GeV extrapolated to the physical point [20] 

0.70 ±0.20 at /i- 1.8 GeV for ~ 0.8 GeV [21] 

0.66 ± 0.07 at fl ^ 2 GeV extrapolated to the physical point [23] 

0.682 ± 0.018 at /i = 2 GeV for m„ = 900 MeV [22]. 

In the CQSM we have 2J'^{0) = 1 for any rriTr at the low scale of the model, since there are only quarks and antiquarks 
to carry the nucleon angular momentum, and they must, of course, carry 100 % of it. Considering evolution one obtains 

2J'3(0) « 0.75 from CQSM at ^ « (1.7 - 2.0) GeV, (32) 

and 27*^(0) « 0.25. The result (32) is — within model accuracy — in good agreement with the lattice data in (31). 

For other aspects in the context of the nucleon spin structure discussed from the CQSM point of view the reader 
is referred to [54]. 



2J'3(0) = < 
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VII. CHIRAL EXTRAPOLATION OF LATTICE DATA 

Let us draw from our study some conclusions which could be of interest in the context of the chiral extrapolation 
of lattice data. Unfortunately, the model cannot provide any insights in this respect concerning the m^-dependence 
of M§ or J^. At the low scale = 2J<3 = 1 which expresses consistently that the total nucleon momentum and 
angular momentum in the model are carried by quarks — irrespective the value of m^r. The model provides, however, 
interesting results for the t-dependence of M^(<), J^{t) and including the normalization of the latter. 

For that let us imagine that for some reason we were able to compute in the CQSM the EMT form factors only for 
m-iT ^ (400 — 500) MeV, and then forced to extrapolate the results down to small TOjt in order to compare to the real 
world — using the model results from the large-m^r region as main guideline. What would we obtain? 

First of all we emphasize that in the CQSM all form factors of the EMT can be well approximated by dipoles. 
Though there is some prejudice in this respect based on the experience with the electric and magnetic form factors of 
the proton, this is an important and non-trivial observation. 

Next, we consider the dipole masses of the form factors M2 (t) and J'^it). From the approximately linear m^- 
behaviour of the respective dipole masses in the rrij^ > (400 — 500) MeV we might have been tempted to assume also 
a linear behaviour in the region m^-region down to the physical point. Such an extrapolation would have resulted in 
a surprizingly accurate estimate for the dipole masses of M^it) and J^{t) at the physical point, see Fig. 8a. 

However, this procedure would have failed in the case of the dipole mass of d'^it). In fact, also this dipole mass 
exhibits an approximately linear m^-behaviour in the region ^ 500 MeV, see Fig. 8a. But presuming that this 
linear behaviour continues to hold also for < 500 MeV would have resulted in a strong overestimate of the dipole 
mass of di{t) at the physical point of about 30%, see Fig. 12a. 

Finally, let us discuss the m,r-dependence of the form factor d^{t) at zero momentum transfer. We observe an 
m,r-behaviour of df = d'^ (0) in the region of m^r > 400 MeV which is linear to a very good approximation. However, 
if we assumed that this linear m^r-dependence continued down to small pion masses we would have underestimated 
the absolute value of df at the physical point by about 25%, and in the chiral limit by 50%, see Fig. 12b. 

Recall that the chiral expansion of d'^ to linear order in — i.e. up to the leading non-analytic contribution (in 
large N^) according to Eq. (26) — approximates the full result rather well up to the physical point. It is an interesting 
result that in our model df exhibits in the region of large m,r > 400 MeV again a linear TO7r-dependence — just as at 
low < 140 MeV, however, with a roughly 30 times smaller slope! 

Of course, we make these observations in an effective chiral theory formulated in the large Nc limit. The situation 
could be different in full QCD simulations performed at finite Nc = 3. However, the model was observed — in spite 
of its caveats like large vs. finite Nc — to provide a qualitatively good description of lattice results on the pion mass 
dependence of the nucleon mass in the region 300 MeV < < 1.5 GeV [37]. Therefore it could be worth to keep in 
mind the lessons we learn here from this model. 




0.2 0.4 0.6 0.8 in GeV 0.2 0.4 0.6 0.8 in GeV 

FIG. 12: (a) The dipole mass M^^p of the form factor (t) vs. m-^ as obtained from the CQSM. The linear niTr -dependence in the region 
> 500 MeV is emphasized, (b) The pion mass dependence of the form factor (t) at zero momentum transfer as obtained from the 
CQSM at low scale. Solid line: The exact result from the CQSM. Dotted line: The chiral expansion of di(m^) to linear order in 
according to Eq. (26). Dashed line: The attempt of a "chiral extrapolation" from the region of large > 400 MeV. In both figures the 
cross marks the physical point. 
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VIII. CONCLUSIONS 



We have presented a study of the nucleon EMT form factors in the large- A^c hmit in the framework of the CQSM 

in the region of pion masses 140 McV < < 1.2 GcV. This work supplements the study of the nucleon EMT form 
factors in the region < < 140 MeV [55] . There is a good reason why these studies have been presented separately. 
The CQSM [38, 39] has been used extensively to study the nucleon at the physical point and in the chiral limit with 
considerable phenomenological success [41-54]. The model owes its success to the facts that it is a QCD-inspired field 
theoretic approach [61-63] which correctly describes chiral symmetry playing an essential role in the description of 
the nucleon. 

Here we studied the CQSM at model parameters corresponding to large pion masses far above the physical point, 
i.e. far away from the chiral limit. Is the model applicable in this situation? First exploratory studies in this direction 
have indicated a positive answer [36, 37]. In order to shed further light to this question we have focussed here on the 
nucleon form factors of the EMT which is a central quantity for any field theoretic approach. 

In our study we have demonstrated that the soliton solutions for large rnvr obtained numerically in [37] correspond 
to stable solitons, which is an important prerequisite for the applicability of the approach. Wc have shown that the 
model results for the energy density, the angular momentum density, and the distributions of strong pressure and 
shear forces are in qualitative agreement with what one expects as the constituents building up the nucleon become 
heavy and the range of the pion cloud diminishes. Namely, the spatial extension of the nucleon shrinks, the energy 
density increases, the absolute values of the forces inside the nucleon increase — resulting in a smaller and more 
tightly bound nucleon. These densities are related to the form factors of the EMT [3] . 

In order to test the model results in a more quantitative way we have compared them to results for the EMT form 
factors from lattice QCD simulations [20-24] performed at lattice parameters corresponding to > 500 MeV. We 
observe a good agreement with the lattice data within an accuracy of (10 — 30)%. Also at the physical point the 
model was observed to work to within a similar accuracy [41-48]. 

The good performance of the model at the physical point and in the region of large pion masses raises the question 
whether the model results could be used as a guideline for the chiral extrapolation of lattice data. Given the generally 
observed accuracy of the model it is clear that it cannot be used as a precision extrapolation tool. However, the model 
results can be helpful in two respects. 

First, our results can be of use as qualitative guidelines for the extrapolation of lattice data. For example, the 
model results indicate that the EMT form factors can be well approximated by dipole fits, and that the dipole masses 
of M^it) and J^{t) exhibit to a good approximation a linear mjr-dependence from = 1.2 GeV down to the 
physical point. But our results indicate also that the chiral extrapolation of the form factor d^{t) from the region 
> (400 - 500) MeV could be a subtle and difficult task. 

Second, the success of the model at large pion masses could provide a justification for applying the "pion cloud" 
idea to the description of nucleon Tip to pion masses of 0(1 GeV) as explored in the finite range regulator approach. 
This approach may provide a useful and precise tool for the chiral extrapolation of lattice data as argued in [33, 34]. 

To conclude, we observe that the CQSM provides a consistent description of the nucleon at large pion masses 
indicating that efi'ective quark and pion degrees of freedom can account for the gross features of the nucleon properties 
also in this regime. Work which may provide further insigths in this direction is in progress. 
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APPENDIX A: ALTERNATIVE DEFINITION OF FORM FACTORS 



By means of the Gordon identity 2Mmu'^°'u = ?i'(«(T"''A„ + 2P")u one can rewrite (1) as (see e.g. in Ref. [5]) 



2 4:M]\[ 



M. 



N 



(Al) 



where AQ'^{t) = M^-^{t), AQ'^{t) + BQ<^{t) = 2 J'^ (t) , C^'^it) = \ d'^''^{t). In this notation the constraints (8) 
read A^{Q) + A'^{Q) = 1 and B*^(0) + B^{Q) = 0. The latter constraint is sometimes rephrased as the "vanishing of 
the total nucleon gravitomagnetic moment" . 
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